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QJ ' Abstract. Colour measurements of quasars are used to constrain the abundance and 

C/^ , properties of intergalactic dust and the related extinction effects on high- 2; sources. For 

I Type la supernovae at z = 1, we derive an upper limit on the possible dimming from 

intergalactic dust of Am = 0.2, ruling out the "grey" dust scenario as being solely 
responsible for the observed faintness of high-z SNIa. 
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1. Introduction 

Scattering and absorption of light by a homogeneous distribution of intergalactic large 
grains giving a nearly achromatic cross-section - so called grey dust - has been proposed 
PI E] as an alternative to the dark energy explanation for the observed faintness of Type 
la supernovae (SNIa) at 2; ~ 0.5 [Sllll 
^ I The presence of a significant component of dark energy has been inferred from 

complementary cosmological tests, e.g., from recent WMAP background radiation 
I measurements [S] in combination with limits on the Hubble parameter from the HST 
^ [ Key Project P and large scale structure measurements from the 2dFGRS survey [Zj. 

However, in Ref. it was demonstrated how the recent CMB anisotropy observations 
can be explained also by an Einstein-de Sitter cosmology with no dark energy if one 
allows for a low Hubble parameter, Hq < 50 km/s/Mpc, the value currently favoured 
by gravitational lensing time delay measurements [Hj- It is further argued that conflict 
with large scale structure measurements can be avoided by suppressing cluster scale 
mass fluctuations assuming a small but non- negligible hot dark matter component. 

Thus, probing the redshift-distance relation through Type la SN observations 
remains essential for establishing and exploring the dark energy of the universe. Ruling 
out alternative explanations to the dark energy interpretation and minimizing systematic 
effects on the measured distances to high-z supernovae is still crucial. In this note, we 
discuss the possibility to constrain the effects of intergalactic grey dust. 
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It has been argued that the amount of intergalactic dust required to make the 
SNIa results compatible with an Einstein-de Sitter universe is disfavored by the far- 
infrared background measured by the DIRBE/FIRAS instruments PO]- However, the 
derived bounds rely on the assumption that optical and UV photons are absorbed and 
re-emitted rather than just scattered by the dust grains. 

It is expected that large dust grains will also cause some differential extinction for 
high-z sources (although less than for smaller dust grains), allowing the possibility to 
distinguish between extinction and a cosmological origin of the dimming. In Ref. it 
was shown that with 1 % relative spectrophotometric accuracy or broadband photometry 
in the wavelength interval 0.7-1.5 fim for high-z SNIa observations, one could reduce 
the uncertainty from grey dust extinction down to Am = 0.02 magnitudes (the target 
error for the proposed SNAP SN survey [T2]). 

In this paper, the differential extinction for high-2; quasars (QSOs) caused by 
extragalactic dust along the line of sight is used to put an upper limit on the possible 
effects on SNIa observations at z = 1 of Am = 0.2 (comparable to current SNIa survey 
errors) for a wide range of models for evolution of dust density and extinction properties. 

2. Extinction by dust at cosmological distances 

For a given emission redshift z^, the attenuation Amd at observed wavelength Aq due 
to dust can be written (see also Ref. [TT]) 



where a(A, Ry) is the dimensionless wavelength-dependent attenuation coefficient where 
we use the parameterization given in [T3] and Dy^z) oc [cr ■ Pdust{z)]~^ is the interaction 
length for scattering or absorption by dust particles in the observed y-band (A = 
0.55 yum). Here, pdust(-2) is the physical dust density at redshift z, a the interaction 
cross-section and the proportionality constant will be related to the mass-distribution 
of the dust particles (since Dv{z) = [cr ■ ndust{z)]~^ where ndnst{z) is the dust number 
density). The cosmology-dependent function h{z) is given by 



Note that we have set c = 1 and h{z) and Dv{z) should be given using the same set of 
units, e.g., Mpc. 

Thus, our parameters are the reddening parameter Ry and the interaction length 
at zero redshift Dov oc (a ■ Pdust)"^- 

The Monte-Carlo simulation program SNOC [TF was used to perform the integral 
in Eq. numerically. Note that the model is valid also for a patchy dust distribution, 
as long as the scale of inhomogeneities is small enough, i.e., 1/\/N <^ 1 where N is the 
number of dust clouds intersected by the light-ray. 

Perfectly grey dust (i.e., with wavelength- independent absorption) can obviously 
mimic the effects of any dark energy component if one allows for arbitrary fine-tuning of 
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the dust distribution. We have investigated whether any reahstic dust model related to 
astrophysical sources are consistent with the available data. We note first that the dust 
density has to be related to astrophysical sources, such as star formation, and second, 
that physically reasonable dust models generally implies a wavelength-dependence in 
the absorption and scattering properties. 

As in Ref. [TT], two dust distributions are considered: pdust oc (1 + z)", where 



Thus, Model A has a constant comoving dust density whereas Model B with 
Pdust(-2^ > 0.5) = Pdust(-2^ = 0.5) implies dust creation at z > 0.5. Models A and B 
are considered to be extreme in opposite directions, i.e. they should capture the full 
range of physical models leading to the production of intergalactic dust. 

In our simulations we explore values of Ry in the interval ranging from 1 to 9 
and the interaction scale-length at zero redshift Dqv in the range 10 — 10 000 Gpc. In 
a [^M = 0.3, = 0.7] -cosmology, z = [0.5,1.0,1.5] correspond to proper distances 
of 1.55, 2.37 and 2.84 Gpc, respectively. Note that since the dust density varies with 
redshift according to Eq. ©, the interaction length will vary correspondingly. 

3. Spectral homogeneity of QSOs 

QSOs have been found to have remarkably homogeneous spectral properties in spite 
of their wide range in brightness. The Sloan Digital Sky Survey collaboration (SDSS, 
[in]) has produced empirical spectral templates based on a couple of thousand objects 
in the redshift range 0.044 < z < 4.789 showing a spectrum-to-spectrum la difference 
of approximately 20 % in the rest system wavelength range 1200 to 8500 A j2ni. HST 
QSO data was used by [23] to extend the spectral template down to 300 A in order 
to obtain full wavelength coverage for all broadband filters (see Sec. E} over the entire 
redshift interval used in this paper. 

4. The QSO sample 

The SDSS and 2dF [T7j Galaxy Redshift Survey will eventually include about 10^ QSOs 
up to redshift z ~ 5. We have used a data sample of 3814 QSOs from the SDSS Early 
Data Release (EDR) JH]- These have redshifts between z = 0.15 and 5.03 distributed 
according to Fig. [H 

The optical magnitudes are measured through five broadband filters (ugriz) with 
limiting point-source magnitudes of 22.3, 22.6, 22.7, 22.4, and 20.5, respectively, for a 
signal-to-noise ratio of 5:1. 

The quoted photometric calibration errors are 3 % in g,r and i, and 5 % in m and z 
PH] • The instrumental response functions are given in Ref. |2ZI where we use the second 
column (respt) corresponding to 1.3 airmasses at APO for a point source. 




3 for all z, Model A 

for z > 0.5 (3 for lower z). Model B 
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Figure 1. Redshift distribution of 3814 QSOs in the Sloan Digital Sky Survey Early 
Data Release. 

We exclude QSOs that are classified as extended to avoid contamination from the 
host galaxy. This cut leaves 3692 of the original 3814 QSOs. Applying the 5:1 signal- 
to-noise ratio limiting point-source magnitude cut leaves 3556 of the point-like QSOs. 
In order to avoid confusion with Seyfert galaxies, we also exclude objects at 2; < 0.5 
[19 , 20 . To avoid putting excessive weight in the analysis on a possible dust distribution 
at very high redshifts (and also to minimize a possible bias from flux limitations) we 
have only considered QSOs with redshift z < 2 in our primary analysis, leaving 2740 
point-like QSOs at 0.5 < z < 2. An alternative method to constrain intergalactic dust 
at 2; > 2 by studying the thermal history of the intergalactic medium is presented in 
Ref. 

5. Method and results 

We divide the QSO sample into redshift bins of size Az = 0.1. In each redshift bin, 
we calculate the mean colour, e.g., (m — 5')obs as a measurement error weighted average. 
Objects that have colours that are more than 2 a off from the mean value are rejected to 
avoid contamination from object with anomalous colours, e.g., QSOs that have suffered 
severe extinction in the host galaxy. 

For each set of parameter values Ry = [1,9], Dqy = [10,5 000] Gpc and redshift, 
we simulate the attenuation due to grey dust and add to the median QSO spectrum 
obtained in Ref. |26j . 
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For a pair of broadband filters X and Y, we define 

A(X - F).,R„Dov = {X- F). -{X- y).,R,,D,, (4) 

where (X — Y)z is tlie observed colour of a QSO at redshift z and {X — y)z,RY,Dov ^^e 
simulated colour for a specific set of dust parameter values, i.e, the colour from the 
template spectrum corrected for differential extinction and convolved with the SDSS 
filter functions. We then calculate the ^ function of Ry and Dqv as 

N 

where the sum is over N redshift bins and V{X — Y)ij is the covariance matrix defined 
as 

V{X-Y),, = a{X)l,+a{Y)l,+6,, HX)l^ + a{Y)l^ + a{X - (6) 

Here, agys denote the systematic errors associated with each broadband filter (see 
Sec. mi, (Ttem the template error (Sec. 16. 2j) and (Jobs the error in the mean observed 
colour in each redshift bin (after the cut of outliers described above). The scatter in 
each redshift bin is approximately 0.2 mag so with the order of 100 QSOs in each redshift 
bin we have cxobs ~ 0.02 mag. 

In Fig. 121 the colour-redshift evolution for a couple of different dust scenarios are 
compared. The upper left panel shows A{u — i) vs redshift for the Model A best fit 
parameters Ry = 9 and Dqv = 100 Gpc and the lower left panel shows A{u — i) 
vs redshift for the Model A best fit parameters Ry = 2 and Dqy = 1000 Gpc. 
These parameters yield x^/dof ~ 0.5, indicating that systematic errors may have been 
overestimated (see Sec. These plots are to be compared with the case of Ry = 4 
and Dqv = 10 Gpc that give high x^- values (i.e., low probabilities) for the same dust 
distribution model (right panels). 

We have also investigated the dimming effect on high-2; sources at different 
wavelengths, in particular the restframe i?-band flux of SNla. In Fig. IHl we have 
combined the results from this investigation with the results from the x^-test. In the 
left panels, the results for A{u — i) assuming the Model A (upper panel) and Model 
B (lower panel) dust distribution are shown. The right panels show the corresponding 
results for A{g — z). The red, dark orange, light orange, and yellow regions indicate 
68, 90, 95 and 99 % confidence level allowed regions from the x^-test. The four contour 
levels labeled [0.2,0.1,0.05,0.02] correspond to the amount of attenuation caused by grey 
dust when integrating over the restframe 5-band of a set of SNIa at a redshift of 2; = 1. 

In Fig. 15 we show combined results from colours A{u — i) and A{g — z). Allowing 
for any value of Ry within the parameter range, we are able to rule out Am > 0.2 
at 99 % confidence level for SNIa at 2; = 1. Assuming intergalactic dust properties 
more similar to galactic dust {Ry ^ 4) we are able to do almost a factor of ten better 
(Am < 0.03). Note also that all our results are consistent with a no dust hypothesis 
since the x^- values flattens out when D^y 00. The small scale features in Fig. Eland 
IHare due to low resolution in the parameter grid. 
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Figure 2. The colour-redshift evolution for different dust scenarios. The upper left 
panel shows A{u — i) vs z for Ry = 9 and Dqv = 100 Gpc (Model A), the lower left 
panel shows A{g - z) vs z for i?y = 2 and Dqv = 1000 Gpc (Model A). The right 
panels show corresponding plots for i?y = 4 and Dqv = 10 Gpc that give very low 
probabilities for the same dust distribution model. 



Redoing the analysis combining A(u — z) and A{g — i) yields very similar results. 

6. Systematic effects 
6.1. Selection effects 

The selection algorithm for the SDSS quasar sample varied throughout the time that 
the EDR observations were obtained, so the objects in this catalog were not found via 
a uniform set of selection criteria Richards et al |22] analysed the limitations on 
the colour range of QSOs and found that for z = 2, the SDSS survey was sensitive 
to reddened QSOs with E{B — V) > 0.54. We conclude that selection effects do not 
significantly affect our ability to compute the average colour of the QSO sample out to 
redshift z = 2. 
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A(u-i), Model A A(g-z), Model A 




Figure 3. Results from restframe B-band attenuation for SNIa at z = 1 and the 
X^-test using A{u — i) (left panels) and A{g ~ z) (right panels). Contour levels labeled 
[0.2,0.1,0.05,0.02] show the amount of attenuation (in magnitudes) while the red, dark 
orange, light orange, and yellow regions indicate 68, 90, 95 and 99% confidence level 
allowed regions from the x^-test. Small scale features are due to low resolution in the 
parameter grid. 

6.2. Spectral evolution of QSOs 

The possible spectral evolution is limited by the nearly Gaussian 20 % variations found 
for the whole sample. As we only consider less than half of the survey's redshift range, 
systematic colour trends with redshift would not exceed ~ 10%. In our analysis, we 
use a systematic uncertainty on the magnitude derived from the median spectrum vs 
redshift of 10 % in m and z and 5% in g,r and i. In general, unless the spectral evolution 
trends exactly cancel the differential extinction from dust, our limits on grey dust are 
conservative in the sense that the reddening that we attribute to possible extinction 
by dust may stem from a redder population of QSOs. Of course, a possible problem 
is the fact that the template spectrum we use is derived from spectra that themselves 
might be affected by dust absorption and extinction. Such an effect would not cause 
any major change in differential evolution of the QSO colours vs redshift. However, 
it should be noted that it could cause a close to constant offset in the colour-redshift 
diagram that might bias the outcome of the x^-analysis. An iterative method where this 
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Figure 4. Combined results from restframe B-band attenuation for SNIa dX z — 1 and 
the x^-test using A(m — i) and /^{g — z). Contour levels labeled [0.2,0.1,0.05,0.02] show 
the amount of attenuation (in magnitudes) while the red, dark orange, light orange, 
and yellow regions indicate 68, 90, 95 and 99 % confidence level allowed regions from 
the x^-test. 



is corrected for by shifting the colour of the template to match the colour of the current 
best fit model at the mean redshift of the QSO sample {z ~ 1.4) yields results similar 
to Fig. 0] with the exception that for Model B, the upper right corner (corresponding 
to high values of Ry and -Dqv) is also excluded. However, the main conclusion that 
Am > 0.2 is excluded for SNIa aX z = 1 remains unchanged. Improved statistics and 
systematic errors in the QSO surveys will improve the control of systematic effects from 
the template spectrum, see Sec. d 



6.3. Host galaxy dust 

Internal reddening may come from depletion by material associated with the torus of 
gas and dust that is thought to surround the central QSO engine [THj. Richards et al. 
[22] estimate that 6 % of the SDSS quasar sample shows signs of reddening by dust, 
probably in the host environment. The computation of average colours for each redshift 
bin takes this into account by an iterative rejection mechanism where QSOs more than 
2 a away from the average colour at that redshift are rejected, see Sec. El Thus, the tails 
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of the distribution are suppressed. Remaining reddened QSOs would tend to make our 
limits more conservative, as the residual host galaxy dust would be counted as possible 
intergalactic dust. 

6.4- Evolution of dust properties 

Although we have assumed a constant value of Ry for the redshift range 0.5 < z < 2.0, 
our derived limits show only moderate dependence on the actual value of the reddening 
parameter. A variation of Ry within the considered range would thus not lead to a 
significant difference in the observed attenuation of high-z SNIa. 

7. Future prospects 

The Sloan Digital Sky Survey has recently made available the Data Release One (DRl) 
with spectroscopic parameters for 17 700 medium- 2; {z < 2.3) and 980 high-z QSOs 
{z > 2.3), including reprocessed data from the Early Data Release. A catalog of DRl 
QSOs is being prepared [24 and we plan to redo the analysis with the full data set 
when the catalog is completed. The increased number of QSOs will of course help in 
decreasing the statistical error. Hopefully, the added data will also be helpful in order 
to improve our understanding of different systematic effects. For example, it should be 
possible to constrain the possible evolution of QSO colours and minimize the systematic 
error from the spectral template by studying the homogeneity of QSO spectra over a 
range of redshifts. 

8. Summary and conclusions 

Studying the redshift evolution of colours of QSOs from the Sloan Digital Sky Survey 
Early Data Release, we constrain the effect from extinction by homogeneous intergalactic 
dust on high-z observations. We compare observed spectra from 2740 QSOs at 
0.5 < z < 2 with simulated spectra obtained by adding the attenuation due to different 
grey dust scenarios to a reference median QSO spectrum. For each dust model, we also 
compute the related extinction effects on high- 2; SNIa. 

For a wide range of models for the dust distribution and extinction properties, we 
are able to rule out Am > 0.2 at 99 % confidence level for SNIa at z = 1. This is 
comparable to current SNIa survey errors and allows us to rule out the "grey" dust 
scenario as an alternative to the cosmological constant explanation for the observed 
faintness of high- 2; SNIa. It should be noted that these results are based on the specific 
parameterization of the wavelength-dependent dust attenuation given in Ref. Jl] and 
that the validity of this specific parameterization for very high values of Ry remains 
untested. 

To constrain the dark energy properties (e.g., the equation of state parameter) with 
future SN surveys such as SNAP, yet better constrains on dimming by grey dust are 
required. Assuming intergalactic dust properties similar to galactic dust (i?y < 4) the 
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present analysis rules out Am > 0.03. However, for a more general dust scenario, larger 
statistical samples and better control of the systematics of the QSO sample are required. 
Hopefully this will be within reach with future data sets from, e.g., the Sloan Digital 
Sky Survey. 
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